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CONCERNING  THE  AUTOMATION  OF  'PROCESSING 
OF  EXPERIMENTAL  DATA  IN 
INVESTIGATING  THE'  IRREGULAR  IONOSPHERE 


|  /Following  is  the  translation  of  an  article  j 

I  :by  V,B*  Gusev  and  T,A.  Gaylit  entitled  "Ob  j 

I  Avtomatlzatsil  Obrabotkl  Eksperimenfc&l snykh  I 

p:rl  Issledovanli  Neregulyarnoy  lono-  1 

sTer-i*  (English  version  above)  In  Yestnlk  I 

.  Mos&ovskogo  Uhiveraltata *  Seriya  III*  FIz*  I 

!  :  Astron,  (Herald  of  the  Moscow  University,  3rd  | 

.Series,  Phys.  Astron,),  Moscow,  I960,  Mo,  1,  | 

pages  39-49*/-  j 

The  article  discusses  the  circuitry  of  an  elec-  I 

,<  ? 

I  fcronic  correlo graph  for  the  calculation  of  the  structure  I 

*  ? 

j  function.,  of  a  random  process  having  a  spectrum  of  0-4  cpa  \ 

I  i 

j  without  previous  recording  on  a  photographic  film  or  mag-  ? 

i  *  f 

netlc  tape.  The  accuracy  of  calculation  is  3$.  j 

|  The  proposed  layout  can  be  used  for  the  ex&minatlbl 

I  of  the  dynamics  of  the  fine  structure  of  the  Ionosphere,  | 
I  and  for  investigation  of  any  random  process  with  a  low-  I 

_ _ _ _ — - - -J 
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INTRODUCTION 


Problems  of  investigation  of  the  fine  structure  of 
the  ionosphere  are  most  fully  solved  at  the  present  time 
by  means  of  correlation  analysis-.  The .  calculation  of  such 


characteristics  of  the  fine ' structure  as  the  dimensions  of  j 
nonuniformitiesj  their  anisotropy,  the  rate  of  their  chaotic 
variations,  direction  and  velocity  of  drift,  involves 
analysis  of  correlation  functions  of  the  reflected  radio 
signal  amplitudes  as  observed  at  three  spatially  separated! 
points  /l/„  I 

.Inasmuch  as  regular  observation  of  the  state  and  ; 
particulafLy  the  changes  of  the  nonuniform  structure  of  j 

the  ionosphere  Is  of  interest,  it  is  necessary  to  automate! 


the  processing  of  the  data  at  least  partially j  the  corre¬ 
lation  function  In  particular  should  be  obtained  automati¬ 
cally,  becaxxse  the  calculation  and  manual  construction  of 
these  functions  is  quite  time  consuming. 

In  the  correlators-described  in  the-  special  litera 
ture  (see,  for  example,  /2/),  the  process  being  Investiga¬ 
ted  is  recorded  on  a  photographic  film  or  magnetic  tape 
and  then  processed  n  times  through  the  correlator  (n  is 
the  number  of  points  of  the  function  to  be  calculated). 

For  solution  of  the  special  problem  arising  in 
j ionospheric  measurements,  a  simpler  device  can  be  con- 
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,  •  V 

sfcructed  which  determines  simultaneously  n  points  of  the 
function  being  calculated  without  a  previous  recording  of 

1 

the  process,  } 

Considerable  simplification  of  the  circuits  can  be 

achieved,  by  using  discrete  selection,. '  j 

,  £ 

In  processings  the  continuous  process  is  replaced  ; 
by  discretely  selected  values  spaced  at  a  time  Interval  t^j 
The  integration  of  the  correlation  function  j 

t  i 

f  x(i)x(t  ±4dt  I 


is  replaced  by  the  summation 


V*  X  {/  tiQ  X  (t  +  -r  t). 


X-0 


A  time  interval  t^  for  which  two  successively  selected  I 

i 

quantities  ac(t)  and  x(tftQ)  do  not  correlate  mutually  is  j 

! 

the  moat  rational  /4,5/*  Let  us  note  that  the  use  of  dls- 

i 

crate  periodic  selection  is  possible  for  random  processes  | 
only,  For  periodic  processes,  such  selection  will,  in  gen¬ 
eral  ,  give  incorrect  results. 


Section  1, 

The  peculiarities  which  arise  in  connection  with 


j construction  of  a  correlation  device  for  study  of  the  fine  j 
| structure  of  the  ionosphere  are:  .  * 

\  1„-  The  part  of  the  frequency  spectrum  of  the  random 

5  | 

|  process  to.be  investigated  is  0-3  ops,  Hie1  radii  of  corre-  | 
|  lation  vary  from  o'ne  to  several  seconds*  -liras  the  range  ] 

i  ■ ,  f 

f from  0.1  to  several  seconds-  is  a  satisfactory  retardation  j 

i  t 

f  time.  | 

I  I 

|  2. ‘The  problem  of  the  averaging  time  of  the  process 

j  3 

|  has  not  been  fully  resolved,.  An  increase  in  the  averaging  j 
j  interval  reduces  the  statistical  error  of  the  function  de~ j 

I  •  ■  I 

|  termination,  but  the  presence  of  large  nonuni form! ties  in  | 

j  the  ionosphere  leads  to  instability  of  the  process,  and  | 

j 

[  plays  an  Important  role  in  the  case  when  the  time  of  iute-| 


gr&tion  is  of  the  same  order  as  the  period  of  a  large  non-j 
|  '  '  * 

I  uniformity.  The  optimum  integration  time  is  three  minutes. j 

|  Here,  the  integration  interval  contains  60-100  correlation! 

j 

I  radii  of  the  process,  the  statistical  error  In  determining! 

I  the  correlation  function  is  10-15#  /3/*  and  the  averaging  j 

!  { 
time  is  a  fraction  of  the  period  of  the  large  nonuniformi-j 

ties  (15-20  minutes),  | 

3.  Since  the  statistical  errors  in  determination  j 

\  of  the  function  may  be  relatively  large  (10-15#)*  the  Pre"| 

|  cision  of  the  apparatus  need  not  be  too  high.  | 
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Section  2.  j 

Any  function  simply  connected  to  the  correlation} 
function  may  be  taken  as  the  function  to  be  calculated*  ; 
Such  a  function,  for  example,  is  the  structure  function,  J 


c  - —  y  ess  2-3*  | 

«••***  ■%  {t\ 

5- » jr—  x"» 

.  .  j«.  —  jc; 


— 2  — i 

X  —X 


|  In  a  similar  manner,  it  is  possible  to  set  up  the  structure 


function  reciprocally t 


nu,u*»  *n> I 

/wM.**  (*  ~~  y<;  } 


i  jxy\-r-  "  »  4 

j  if  %  y*  s*  "  ®  2'r|!—- 

j  If  x  and  2_ are  amplitudes  of  a  reflected  signal  observed 

| 

j  at  two  separate  points  in  space,  the  assumption  of  equality 
of  the  averages-  and  dispersions  of  these  quantities  corre¬ 
sponds  to  the  assumption  of  a  steady  and  uniform  statisti-l 
cal  ensemble  of*  measured  amplitudes.  For  this  condition,  j 
;  which  is  usually  assumed  to  hold,  the  norms  of  the  8 true - 
|  ture  functions  of  the  amplitudes  are  equal  when  they  are 
I  observed  at  points  separated  in  space.  Therefore,  if  three 


receiving  channels  with  equal -amplification  at  the  operat- 
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k 


i 


£ 

Ing  frequency  are  used,  there  is  no  need  for  normalization f 

I 

of  the  structure  functions  for  calculating  all  average  ! 
parameters  of  the  motion  and  the  geometry  of  the  diffraction 
pattern,  except  for  the  correlation  radius.  If  necessary,  \ 

i 

the  function  may  be  normalized  with  the  aid  of  the  corre-  ; 

lator  to  be  described  (see  -Section  5).  ; 

The  structure  function  is  better  suited  for  auto-  f 

1 

matic  computation  than  the  correlation  function  for  the  | 

following  reasons :  1 

$ 

•i 

1.  Fewer  and  simpler  operations  are  required  for  J 

calculation  of  the  structure  function.  | 

| 

2.  The  calculation  of  the  correlation  functions  j 

I 

requires  very  high  accuracy  in  some  operations.  Indeed,  j 

1 

the  difference  xy_~  xy  for  a  random  process  is  consider-  f 
ably  smaller  than  each  of  the  above  terms  xy^  and  xy .  The  j 


j  storage,  multiplication,  integration,  and  subtraction  er-  f 

1  _  __  | 

|  rors  apply  precisely  to  the  quantities  xyr  and  xy.  The  j 

s  I 

|  relative  error  of  the  difference  may  be  several  times  { 

|  greater  than  the  errors  of  the  subtrahend  and  minuend .  In  j 

I  calculating  the  structure  function,  only  storage  and  sub-  \ 

I  ■  f 

!  traction  errors  can  give  this  effect,  and  higher  require-  * 

j  ■  ~  ■  j 

\  merits  are  set  up  for  the  precision  of  these  operations.  j 

1  | 

|  We  can  also  use  a  simpler  function  —  the  average  \ 

I  \ 

f  modulus  of  the  difference  jx  —  xp„  The  advantages  of  this  | 
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| function  are  obvious.  A  disadvantage ,  however,  is  that  this! 

i; function  is  in  known  relationship  to  the  correlation  func-  j 

tion  only  for  the  normal  and  Rayleigh  distributions.  > 

«  ’  * 

\  Below,  a  device  for  setting  up  the  structure  func-  j 

| tion  of  a  random  process  is 'discussed,  but  its  principles  j 

I  I 

are  also  applicable  to  computation  of  the  correlation  fune~- 

!  '  '  *  *  .  \ 

\  tion*  '  ] 


j  Section  3*  ■  .  * 

I 

The  process  of  computing  a  structure  function  may  | 

.  -j 

|  be  divided  into  the  following  steps*  time  delay  of  the  .  | 

|  process,  subtraction,  squaring.  Integration  and  display.  I 

Figure  1  shows  the  block  diagram  of  the  equipment, j 
-and  Fig,  2  shows  the  shape  of  the  voltages  at  different  | 
points  of  the  ‘circuit,  J 

1.  A  capacitor  is  used  as  a  storage  device.  Switching i 
of  the  capacitor  and  all  other • commutations  are  effected • 

|  using  a  type  Shi -26  stepping  switch.  The  capacitor,  which  | 

1  is  connected  to  the  moving  contact  of  this  selector,  is  j 

f  n 

1  charged  through  the  first  stationary  contact  to  the  voltagt 
I  | 

|  x(t)  being  examined  at  the  start  of ' the  cycle*  Then  the  j 

*  i 

j  capacitor  is  connected  in  series  with  the  input  of  the  j 

|  subtracting  unit  at  the  moments  of  time  t+t^,  t+Tg,  etc-.,  J 
|  thereby  supplying  constant-amplitude  pulses  till  the  be-  J 
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l ginning  of  the  next  cycle.  In  the  next  cycle,  the  constant  j 
value  to  he  stored  is  approximately  equal  to  the  voltage  j 
j being  studied  at  the  moment  of  time  tttQ,  etc. 

i  .  F 

|  Tine  time-variable  voltage  being,  studied  is  applied  j 

|  to  the  second  input  of  the  subtracting  unit,  again  by  means 
'  ; 

I  .of  the  stepping'  switch.  .  ; 

1  ■  > 

|  lfhus,  the  input  of  the  subtracting  device  receives 

I  '  f 

}  two  voltages  with  different  time  shifts.  Therefore,  n  • 

|  •  ~~  f 

I  points  of  the  correlation  function  can  be  found  by  one  sub-r 
tracting,  and  squaring  unit,  1 


7X  “*■ 


* 
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j  Pig,  1.  .Block  diagram  of  an  automatic  structure-function  J 
1  computer*  1}  storage  elements  2)  subtracting  and  squaring  j 
!  unit,;  3)  integrators  i  4)  direct-current  amplifier  j  5)  in-  j 
I  dicator  tube  with  photo  attachment.'  \ 


The  square  of  the  unknown  voltage  is  supplied  to  j 

£ 

the  integrators  through  another  bank  of  the  stepping  switcl 


If 
in  such  a  manner  that  the  results  of  subtraction  and  squ&r-  \ , 

|  ing  which  belong  to  one  value  of  the  time  shift  are  added  j 

in  each  integration.  After  completing  the  integration,  the  ;. 

|  voltage  is  taken  from  the  integrators  by  the  stepping  swit<jh 

and  supplied  via  .the  direct-burrent  amplifier  to  the  verti-  j 

i 

cal  deflecting  plates  of  the' indicating  oscilloscope.  Here^ 

♦  •.  •  'I 

a.  discrete  sweep  is  supplied  to  the  horizontal  deflection  } 

! 

|  plate;  Thus  the  structure  function  is  displayed  on  the  tub 
|  screen,  and  can  be  observed  visually  or  photographed. 

Let  us  now  consider  the  details  of  the  separate 

’  .  ‘  ,  .S'"'  g 

blocks* 

1.  The  storage  errors  result  from  incomplete  charg 

$ 

ing  of  the  capacitor  and  partial  discharge  on  connection  j 


to  the  subtracting  unit. 

The  conditions  for  the  minimum  storage  errors  are* 


«!/<£  R. 


i**jp 

l)  charging!  2)  discharging.  | 

where  At  Is  the  duration  of  the  stepping-switch’s  contact  { 

time  and  depends  on  the  velocity  of  its  motion!  n  is  the  j 

number  of  points  of  the  function  to  be  calculated,  j 

The  maximum  frequency  of  the  process-  which,  is  not 

distorted  by  the  selector  and  the  duration  of  the  delays 
whl ch  can  be  achieved  with  the  apparatus  are  dependent  .uni 
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/•  <* 

t 


the  velocity  of  movement.  ' 

In  our  case*  Csfcor  »  25^F,  Kdig  9*5  megohms,  and 
Rch  as  500  ohms.  These  parameters  ensure  an  accuracy  of  0.5 
for  a  stepping-switch  rate  of  10-20  steps  per  second.  The 

j upper  passed  frequency  Is  2-3  on s,  the  minimum  delay  is 

t 

|  0,05  seconds,  and  the  maximum  is  3  seconds. 

2,  We  obtain  the  absolute  value  of  the  differ¬ 
ence  at  the  output  of  the  subtractor,  because  the  dynamic 
range  doubles  when  the  modulus  is  fed  to  the  input  of  the 
square-law  generator.  The  error  of  the  subtractor,  referre| 
to  the  maximum  value  of  the  difference.  Is  0.5$. 

3,  Squaring  is  carried  out  on  alternating  current 
Triangular  pulses  with  period  and  amplitude  proportional 
to  the  input  'voltage  are  shaped  In  the  square-law  genera t of . 
The  error  of  this  generator,  referred  to  the  whole  range 
of  its  operation,  Is  1$. 

•  The  circuits  of  the  subtracting  and  squaring  de¬ 
vices  are  shown  In  Fig.  3* 

4,  Rather  large  RC  networks  (RC  ~  104  sec)  are  use| 
as  Integrators  in  the  circuit.  The  use  of  such  simple  in- 

1  tegrators  is  possible  because  due  to  the  switching  of  the 
selector  switch,  the  integration  time  is  a  fraction  of  the 
time  during  which  the  process  is  being  examined. 

If  the  process-examination  time  Is  2-3  minutes. 
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the  operation  time  of  the  Integrator  Is  only  about  5-7  sec-j 

2 

onds.  For  RC  **■  104  sec,  the  error  of  Integration  Is  2.5-3$.j 
Tube  circuits  should  be  used  to  obtain  high  precision 
5,  As  noted  above ,  the  indicator  is  an  oscilloscope?,, 

|  and  the  calculated  function  is  displayed  as  a  series  of  | 


v 

■discrete  points  on- its  screen 


6 .  The  possibility  of  error  due  to  incorrect  trans-j 

{ 

mission  of  the  voltages  by  the  stepping  switch,  especially  I 


for  values  v,  should  be  noted.  In  order  to  eliminate 


errors  of  this  kind,  the  condition  of  the  selector  contact^ 
must  be  checked. 

Section  4. 

A  noise  generator  was  used  for  a  final  check  of 
the  operation  of  the  apparatus j  this  generated  noise  with 
a  known  correlation  function  which,  in  its  spectral  corapo- 

[ 

sition  and  distribution,  approximated  the  random  variation! 

of  the  amplitude  of  a  signal  reflected  by  the  ionosphere,  j 

The  apparatus  was  used  for  computing  the  structural  fune-  j 

1 

tlon  and  the  average  difference  modulus  of  the  random  pro¬ 
cess,  The  noise  amplitude  was  simultaneously  recorded  on 
photographic  film  and  tabulated,  and  the-  same  function  eomf 

i 

puted  on  a  type  SASL  computer.  I 

In  order*  to  separate  the  instrumental  error  from  f 

_ I 
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Pig,  4,  Structure  function  of  a 
random  process;  0  —points  ob¬ 
tained  with  the  correlator;  +*  — 
calculated  points. 


Fig,  5*  Me  an  mo  dul us  of  d i f fe r- 
ence  of  a  random  process;  0  — 

■''v 

points  obtained  with  the  corre¬ 
lator;  Hr  —  calculated  points 
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the  statistical  error  due  to  the  finite  discrete  selection., 


v&r.  jnwMcr- 


discrete  values  of  the  amplitude  selected  by  the  apparatus! 

were  brlghtend  for  ,  photography  of  the  noise  ;  then  dis~f 

* 

crete  amplitude  values  selected  by  the  instrument  were  read 
Figure  4  shows  structure  functions  of  a  random  prol- 
cess  obtained  with  the  apparatus  and  computed  with  an  SASLjj, 
Figure  3 ' shows  the  corresponding  mean  moduli  of  the  dif-  j 
f erence .  Such  a  comparison  was  repeated  several  times.  The| 
instrumental  error  amounts  to  approximately  2$  for  the  1 

mean  modulus,  and  3-4^  for  the  structure  function,  when  j 
It  is  referred  to  the  maximum  calculated  value.  It  should  1 
be  noted,  that  the  tabulation  error  Is  included  in  this  er-l 

ror.  I 

Figure  6  shows  several  structure  functions  of  thdj 
amplitude  of  a  signal  reflected  from  the  ionosphere  as  ob-j 
tained  with  the  correlator.  I 


Section  3» 

•  Owing  to  the  application  of  discrete  selection,  i 
is  possible  to  obtain  within  the  same  time  Interval  four 
structural  functions  for  three  signals  with  one  function 
generator.  This  is  necessary  in  investigations  of  the 
ionosphere. 

This  is  achieved  by  alternately  feeding  signal 


I 


Famp: 


amplitudes  from  all  three  observation  points  to  the  input 

of  block  2  (Fig.  1)  with  the  aid  of  the  selector  switch  in 

such  a  manner  that  three  amplitudes  with  all  time  shifts 

are  on  the,  same  time  axis*  The  separation  of. these  values 

for  the  respective'  integrators  is  accomplished  by  the  stepj* 

( 

ping  switch.  ■  I 

?  :  j 

The  norm  of  the  structure  function  2d"  is  equal  toj 

(x  ~  xtY'  for  T  for  which  E (t)  »  0.  The  delay  of  the  pro-! 

cess  by  a  time  T  4~5  times  greater  than  the  radius  of  corJL 

I 


relation  is  achieved  by  an  additional  element  consisting 
|  of  storage  capacitors  and  a  stepping  switch  which  has  a  1 
I  rate  1/4  to  1/J>  that  of  the  other  selectors*  Subtraction 
f  £'«nd  squaring  of  this  term  is  obtained  by  using  the  blocks'  } 

s 

described  above. 

In  this  manner*  one  of  the  points  on  the  indicate 
(for  example*  the.  last  one)  vj3.11  correspond  to  2 o'"  or  the 
unit  level  of  the  function  being  computed, 

CONCLUSIONS 

|  The  operating  principle  of  the  automatic  eaulp-  < 

j  . 

|  ntent  discussed  above  is  applicable  'to  the  processing  of 
random  process  data  with  a,  frequency  spectrum  of  0-3  cps,  f 
The  spectrum  can  be  broadened  by  increasing  the  speed  off 
I  the  stepping,  switch  and.  by  appropriate  variation  of  the 
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(equipment  parameters.  The  error  of  the  equipment  is  2-4$.  j 

The  electrical  circuitry  of  the  system. is  simple j  f 

( 

jit  contains  eleven  electron  tubes,  not  counting  the  indicator 
[tube  and  the  power  supply.  ’  j 

There  are  no  special'  mechanical  devices  requiring  I 
[precision  fabrication,  such  -as  tape  transports.  * 

'  ■  t 

The  computing  time  of  the  function  coincides  with  j 
the  time  of  observation  of  the  process ;  there  is  no  time  j 

loss  for-  preliminary  fixation;,  therfore  the  variations  of  1 

■  ■  '  '  1 
the  .-correlation  characteristics  of  the  reflected  field  can  j 

fbe  examined.  '  j 
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